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Abstract 



We numerically investigate the effects of inhomogeneities in the energy 
spectrum of aperiodic semiconductor superlattices, focusing our attention on 
Thue-Morse and Fibonacci sequences. In the absence of disorder, the corre- 
sponding electronic spectra are self-similar. The presence of certain degree 
of randomness, due to imperfections occurring during the growth processes, 
gives rise to a progressive loss of quantum coherence, smearing out the finer 
details of the energy spectra predicted for perfect aperiodic superlattices and 
spurring the onset of electron localization. However, depending on the de- 
gree of disorder introduced, a critical size for the system exists, below which 
peculiar transport properties, related to the pre- fractal nature of the energy 
spectrum, may be measured. 
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I. INTRODUCTION 



One of the most appeahng motivations for the experimental study of aperiodic super- 
lattices (ASLs) arranged according to Fibonacci [0] and Thue-Morse ||^ sequences is the 
theoretical prediction that these systems exhibit a highly-fragmented energy spectrum dis- 
playing self-similar patterns From a strict mathematical perspective, it has been 
proven that the spectra of both Fibonacci and Thue-Morse lattices are Cantor sets in the 
thermodynamic limit [^,0. From an experimental point of view, however, two major limi- 
tations appear to validate their peculiar fractal nature. 

In the first place, it is not possible to fabricate perfect aperiodic structures. Although 
X-ray diffraction studies show that the characteristic structural order of ASLs is preserved 
under moderately large growth fluctuations |^. the way this robust aperiodic order can 
influence the transport properties of actual, defective ASLs has not been received a proper 
treatment in the literature yet. In this sense, the observation of inhibition of vertical trans- 
port in periodic superlattices with intentional disorder in agreement with the theory of 
localization in one-dimensional disordered systems, suggests the possible existence of a com- 
petition between the long-range aperiodic order and the unintentional short-range disorder 
and opens new perspectives |I0| . 

In the second place, even in the most favorable experimental conditions, only finite 
arrangements with a limited number of layers can be manufactured. In this sense, the 
observation of fragmentation patterns in the energy spectra of short ASLs using different 
experimental techniques [0, encourages further analyses as to whether fractal-like spectra 
are to be expected in actual ASLs with an increasingly large number of layers. 

Two important questions then follow quite naturally: First, what are the effects of 
unintentional disorder in the splitting scheme of the energy spectrum of ASLs? Second, 
what are the finite size effects on their fractal-like properties? In this paper we address 
these questions by means of the study of the phase diagram, the locahzation length and 
the bandwidth-scaling of the energy spectra. Our results indicate that small fluctuations in 
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the sequential deposition of layers considerably smear out the self-similarity of the energy 
spectra on increasing the system size. Thus, we conclude that fractal-like spectra with a 
richness of finer details such those found by numerical analyses of perfect ASLs are not to 
be expected in large ASLs, but quantum coherence is strong enough in short systems to 
give rise to a measurable hierarchical set of subminibands in the electronic spectra and to 
influence its related transport properties. 



II. MODEL 

We consider quantum-well based GaAs-Gai_^Al^.As superlattices with the same barrier 
thickness b in the whole sample. The height of the barrier for electrons is given by the 
conduction-band offset at the interfaces. We take the origin of electron energies at the GaAs 
conduction-band edge. The thickness of each quantum-well is Azn — b = — — b, 
being the position of the center of the nth barrier. We will focus on electronic states close to 
the bandgap and neglect nonparabolicity effects hereafter, so that a one-band Hamiltonian 
suffices to describe those states. Physical magnitudes of interest can be easily computed 



using a transfer-matrix formalism in this simple picture ||12[- In particular, we can obtain 
subband energies under periodic boundary conditions and the Lyapunov coefficient, which 
measures the inverse of the localization length in units of the superlattice period. We will 
consider Thue-Morse superlattices (TMSLs) and Fibonacci superlattices (FSLs) since they 
have been already constructed starting from two basic building blocks A and A by means of 
molecular beam epitaxy ||l],|2| . In our model, we take A (A') consisting of a quantum- well of 
thickness a (a') and a barrier of thickness b. ASLs are generated using the following inflation 
rules: A AA', A' A' A for the TMSL and A AA', A' A for the FSL. In this way, 
finite and self-similar ASLs are obtained by n successive applications of these rules, with 
N = 2^ wells in the TMSL case and = F„ wells in the FSL case. The Fibonacci numbers 
are generated from the recurrence law F„ = -F„_i -|- -F„_2, starting with Fo = Fi = 1. 

Interface roughness appears during growth in actual ASLs: Protrusions of one semicon- 
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ductor into the other cause in-plane disorder and break translational invariance parallel to 
the layers. Because the in-plane average size of these defects depends on the growth con- 
ditions and it is unknown in most cases, one is forced to develop a simple approach. We 
carry out such an approach describing local excess or defect of monolayers by allowing A2;„ 
to fluctuate uniformly around the nominal values a + b or a' + b. This approach enables us 
to restore the translational symmetry parallel to the layers, thus facilitating computations. 
Our approximation should be valid whenever the mean-free-path of electrons is much smaller 
than the in-plane average size of protrusions as electrons only see micro-quantum- wells with 
small area and uniform thickness. Therefore, each micro-quantum-well presents a slightly 
different value of its thickness and, as a consequence, resonant coupling between electronic 
states of neighbouring GaAs layers is decreased. To get an accurate description of electron 
dynamics, average over all possible configurations of disorder is indeed required because the 
number of interface defects as well as their mean thickness vary from layer to layer. For 
definiteness we take Azn — a{l + Wen) + b or Azn — a'{l + Wen) + where is a positive 
parameter measuring the maximum fluctuation and e^'s are distributed according to a uni- 
form probability distribution P(e„) = 1 if |e„| < 1/2 and zero otherwise. Note that e„ is a 
random uncorrelated variable, even when the lattice is constructed with the constraint that 
the mean values of A^;^ follow the aperiodic sequences. 



III. RESULTS AND DISCUSSIONS 

We have studied GaAs-Gao.65Alo.35As ASLs with a 6 = 32 A and a! — 35 A. In this 
case the conduction-band offset is 250 meV and the effective masses are m^^^ — 0.067m 
and mQa,AiAs — 0.096m, m being the free-electron mass. In a periodic superlattice with 
a' = a = 32 A only one allowed miniband hes below the barrier. We restrict ourselves to the 
study of the fragmentation of this miniband when aperiodicity is introduced. Averages over 
possible configurations of disorder comprised a number of realizations varying from 50 up to 
100 to test the convergence of the computed mean values, and this convergence was always 
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satisfactory. In order to compare with possible experimental situations, we have considered 
W ranging from up to a maximum of 0.08. This value amounts to having maximum 
protrusion thicknesses of half a monolayer on average. 

Figure |l|(a) shows the dependence of the energy spectrum structure on the amount of 
disorder, W, present in a FSL with N = Fn = 144 wells. The energy is measured from 
the conduction-band edge in GaAs. For a perfect (W = 0) FSL the overall structure of the 
energy spectrum is characterized by the presence of four main subbands. Inside each main 
subband the fragmentation pattern follows a trifurcation scheme in which each subband 



further splits from one to three subsubbands ITSlJT^ . Therefore, the energy spectrum of 
perfect and finite FSLs presents distinct pre-fractal signatures. The situation changes when 
randomness is introduced. In fact, although the tetrafurcation pattern of the perfect FSL 
still remain, the finer details corresponding to successive steps in the hierarchical splitting 
scheme are progressively smeared out on increasing the disorder due to growth fluctuations. 
Figure |I|(b) shows the dependence of the spectrum structure with the degree of disorder 
present in a TMSL with = 128 wells. For a perfect system the fragmentation scheme 
agrees with that previously discussed by Ryu et al. [Q in the framework of the Kronig- Penney 
model, and displays pre-fractal signatures as well. The effects due to unintentional disorder 
are similar to those shown in Fig. ||(a) for a FSL. In a previous paper we have demonstrated 
that pre-fractal signatures of FSLs can be properly described by considering the resonant 



coupling between electronic states of nearest-neighbor building blocks |TJ]. In this sense, 
the above results reinforce this conclusion for they show that the presence of the short-range 
disorder reduces the resonant coupling between quantum-wells and, consequently, weaken 
the physical mechanism giving rise to the self-similar pattern. 

To investigate the effects due to the competition between long-range aperiodic order and 
short-range disorder on the transport properties, we have evaluated the localization length, 
£, as a function of the energy, a magnitude which can be readily determined numerically 
within the transfer matriz formalism [|T2|,|T3[. Figure ^(a) compares the localization length 



for a FSL with = Fu = 144 wells in two different cases. For the perfect {W = 0) 



case we obtain a very spiky structure. Each peak corresponds to a quasi-level, and most of 
them completely extends through the superlattice, its localization length being one order of 
magnitude greater than the system length. Note that the distribution of peaks reflects the 
overall fragmentation of the energy spectrum. On the contrary, when growth fluctuations 
are introduced, the localization length distribution becomes smoother and its value remains 
always smaller than the superlattice length, clearly revealing the onset of localization effects. 

In Fig. 1(b) the same comparison is made for a TMSL with N = 128 wells. Although the 
gross features of the plot, and their corresponding physical interpretation, are completely 
analogous to the case of the FSL, two interesting remarks are in order. In the first place, 
we observe that, in the perfect TMSL, there exists a significant number of states whose 
localization length is several orders of magnitude greater than the system size. This fact 
is related to the presence of latticelike states in the TMSL 0. The occurrence of such 
states does not takes place in the FSL. In the second place, the mean value of £ is lower 
for the defective TMSL, hence indicating that the effects of localization are more intense 
in FSLs than they are for TMSLs with the same amount of disorder. A topic which has 
deserved some attention recently concerns the comparison between the transport properties 
in different kinds of aperiodically ordered structures. Our results indicate that, in presence 
of the same degree of growth fluctuations, TMSLs should exhibit better transport properties 
than FSLs. Finally, notice that the localization onset is more pronounced at the edges of 
the energy spectrum, meanwhile in the central regions of the spectrum (about 110 meV) the 
localization length almost equals the system size. 

From Figs. ^ and ^ we can state that the general effect of disorder in ASLs is twofold. 
Firstly, it masks typical pre-fractal features smearing out the hierarchical splitting scheme 
taking place in perfect ASLs. This effect is independent of the precise nature of the under- 
lying aperiodic order, this is to say, quasiperiodic (FSL) or self-similar but not quasiperiodic 
(TMSL). Secondly, the overall decrease of the localization length of the electronic states 
indicates that localization effects induced by random fluctuations are starting on. 

To get a more complete understanding of the competition between long-range aperiodic 
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order and random fluctuations, we have performed a bandwidth-scaling study [jT6[. To 



this end, we define the equivalent bandwidth 5* as the sum of all allowed energy regions, 
thus being nothing but the Lebesgue measure of the energy spectrum, and compute it as a 
function of the number of wells A^. The obtained results are presented in Fig. ^. Let us start 
discussing the scaling of the Fibonacci system. The equivalent bandwidth in the perfect 
system [Fig. §(a)] decays as a power law of the form S{N) = SqN"^, where 5*0 ~ 40meV 
is close to the bandwidth of the periodic superlattice and P ~ 1/4, within our numerical 
accuracy. According to earlier works , such a behaviour is characteristic of a fractal-like 
energy spectra, becoming a Cantor set in the limit oo. However, in the presence of 

moderate fluctuations {W < 0.05), the equivalent bandwidth decays faster than a power law, 
as seen in Fig. ^(b) for W = 0.05. The deviation from a power law is a natural consequence of 
the breaking of the self-similarity of the electronic spectrum |16|. Using a least square fit we 



have found that 5'(A^) = SoN~^ exp[—a{W)N], where a(W^) ~ 2W'^. From this expression 
two consequences can be drawn. First, the presence of the exponential factor indicates that, 
when disorder is introduced in an aperiodic system, its spectrum is not longer described 
in terms of a pure singular continuous component, whose presence is characterized by the 
presence of the power-law factor N~^. Second, the competition between aperiodic long- 
range order and local disorder can be properly estimated from the deviation of S from a 
pure power-law behaviour. Thus, localization effects dominate the electron dynamics in 
ASLs whose length is larger than the threshold length given by Nth{W) ~ l/a{W). This 
value is about N^h — 200 for W = 0.05. Interestingly the deviation from power law decay 
becomes to be appreciable for moderately long systems, so that short imperfect ASLs will 
show energy spectra quite similar to that corresponding to a perfect ASL. This explains the 
success of previous experimental works to detect self-similar features in actual ASLs pT| . 

Let us now consider bandwidth scaling in TMSLs. As it was formerly discussed by 
Riklund et al. the electronic properties of Thue-Morse lattices lie between those of 
quasiperiodic (Fibonacci) and usual periodic ones. This point has been further elaborated 
by Ryu et al. [Q who have shown that, attending to its scaling behaviour, the spectrum 



of TMSL can be decoupled into two components. One component scales as N~" with 
/? — > 0+, and corresponds to wave functions with a marked degree of extension (latticelike 
states). The other one also scales according to a power-law but /3 > in this case, and 
the corresponding charge density distribution of the related wave functions present large 
spatial fluctuations. According to this, the scaling behaviour shown in Fig. |^(a) can be 
easily understood. For short TMSLs the contribution of both kinds of critical states to the 
total bandwidth is essentially the same, for both spread uniformly over the whole sample in 
an analogous manner. But as the number of wells progressively increases, the contribution 
to S due to the second kind of wave functions decreases as a power law, whereas latticelike 
states give an almost constant contribution to S. As a consequence, only these states 
contribute significantly to the value of S in the thermodynamic limit. The maximum number 
of wells considered in our computations has been = 4096, since the spectrum becomes 
too fragmented to compute accurately S for larger systems. Conversely, the equivalent 
bandwidth S decays exponentially when randomness takes place, as shown in Fig. |^(b). 
It is worth mentioning that asymptotic values of S are almost the same for both kinds of 
disordered ASLs, indicating that the particular kind of long-range order of the underlying 
structure is immaterial with regard to their spatial extend. 



IV. CONCLUSIONS 

In summary, we have proposed a realistic model to study actual aperiodic superlattices 
by allowing the quantum-well thicknesses to fluctuate around its nominal values in order to 
take into account interface roughness. Our results indicate that moderate fluctuations in 
the sequential deposition of layers have significant effects on both the energy spectrum and 
the spatial extension of wave functions. The fractal-like nature of an arbitrary spectrum 
is determined by two complementary features. In the first place, the energy spectrum 
becomes more and more fragmented as the ASL length grows. The physical origin for this 
fragmentation stems from resonant tunneling effects between electronic states of neighboring 
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quantum- wells (short-range effects) . In the second place, the sphtting scheme of the energy 
spectrum must display a self-similar pattern. The physical origin for this self-similarity can 
be traced back to the structural self-similarity of the superlattice itself which, in turn, is 
imposed by the aperiodic ordering of the system (long-range effects) . Taking both facts into 
account we conclude that the purported robustness of the aperiodic order present in ASLs 
does not suffices, by its own, to guarantee the fractal-like nature of the energy spectrum in the 
presence of disorder, because of the main effect of growth fluctuations is precisely to weaken 
the resonant coupling between electronic states. Hence, albeit the structural aperiodic order 
is preserved in the presence of moderate fluctuations, the self-similarity related to it can not 
be properly expressed, in its finer details, in the energy spectrum due to the loss of quantum 
coherence as a consequence of short-range effects. 

The relative importance that this competition between long-range order and short-range 
disorder has on the transport properties depends critically on the length of the system. The 
values of the localization length we have obtained indicate that wave functions no longer 
spread over the whole ASL, as they do in the perfect case, but their degree of extension 
amounts a significant fraction of the system size in contrast with the usual view of localized 
states extending just over a few wells. Thus, as a final conclusion, we can state that distinc- 
tive features of a fractal-like spectrum can be experimentally observed in ALSs of practical 
interest, whose length is smaller than a threshold length. The value of this threshold length 
depends on the heterostructure quality attained during the growth process. 
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FIGURES 

FIG. 1. Phase diagram for a (a) GaAs-Gao.65Alo.35 As FSL with N = Fu = 144 wells and (b) 
GaAs-Gao.65Alo.35As TMSL with = 128 wells. Averages were taken over 10 realizations of the 
super lattice. 

FIG. 2. Localization length for (a) a perfect {W = 0) GaAs-Gao.65Alo.35 As FSL with 
N = Fii = 144 wells (solid line) as compared with a disordered FSL {W = 0.05) of the same 
length (dashed line), and (b) the same comparison for a GaAs-Gao.65Alo.35 As TMSL with A'' = 128 
wells. Averages were taken over 50 realizations of the superlattice. 

FIG. 3. Equivalent bandwidth 5 as a function of the number of wells in (a) perfect W = and 
(b) imperfect W = 0.05 GaAs-Gao.65Alo.35 As FSL (solid lines) and TMSL (dashed lines). Averages 
were taken over 50 realizations of the superlattice. 
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